home *** CD-ROM | disk | FTP | other *** search
/ Cream of the Crop 1 / Cream of the Crop 1.iso / PROGRAM / NRPAS13.ARJ / TUTEST.DEM < prev    next >
Text File  |  1991-04-29  |  1KB  |  58 lines
  1. PROGRAM d13r5(input,output);
  2. (* driver for routine TUTEST *)
  3. CONST
  4.    npts=3000;   (* npts must be > mpts for TUTEST *)
  5.    mpts=600;
  6.    eps=0.03;
  7.    var1=1.0;
  8.    var2=4.0;
  9.    nshft=10;
  10. TYPE
  11.    glnarray = ARRAY [1..npts] OF real;
  12.    glmarray = ARRAY [1..mpts] OF real;
  13.    glnparray = glnarray;
  14. VAR
  15.    glinext,glinextp : integer;
  16.    glma : ARRAY [1..55] OF real;
  17.    gliset : integer;
  18.    glgset : real;
  19.    fctr1,fctr2,prob,t : real;
  20.    i,idum,j : integer;
  21.    data1 : glnarray;
  22.    data2 : glmarray;
  23.  
  24. (*$I MODFILE.PAS *)
  25. (*$I AVEVAR.PAS *)
  26.  
  27. (*$I GAMMLN.PAS *)
  28.  
  29. (*$I BETACF.PAS *)
  30.  
  31. (*$I BETAI.PAS *)
  32.  
  33. (*$I RAN3.PAS *)
  34.  
  35. (*$I GASDEV.PAS *)
  36.  
  37. (*$I TUTEST.PAS *)
  38.  
  39. BEGIN
  40. (* generate two gaussian distributions of different variance *)
  41.    gliset := 0;
  42.    idum := -1773;
  43.    fctr1 := sqrt(var1);
  44.    FOR i := 1 to npts DO data1[i] := fctr1*gasdev(idum);
  45.    fctr2 := sqrt(var2);
  46.    FOR i := 1 to mpts DO data2[i] := (nshft DIV 2)*eps+fctr2*gasdev(idum);
  47.    writeln;
  48.    writeln('Distribution #1 : variance = ',var1:6:2);
  49.    writeln('Distribution #2 : variance = ',var2:6:2);
  50.    writeln;
  51.    writeln('shift':7,'t':8,'probability':16);
  52.    FOR i := 1 to (nshft+1) DO BEGIN
  53.       tutest(data1,npts,data2,mpts,t,prob);
  54.       writeln((i-1)*eps:6:2,t:10:2,prob:11:2);
  55.       FOR j := 1 to npts DO data1[j] := data1[j]+eps
  56.    END
  57. END.
  58.